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Abstract 

Background: Inherited retinal disorders are clinically and genetically heterogeneous with more than 150 gene 
defects accounting for the diversity of disease phenotypes. So far, mutation detection was mainly performed by 
APEX technology and direct Sanger sequencing of known genes. However, these methods are time consuming, 
expensive and unable to provide a result if the patient carries a new gene mutation. In addition, multiplicity of 
phenotypes associated with the same gene defect may be overlooked. 

Methods: To overcome these challenges, we designed an exon sequencing array to target 254 known and 
candidate genes using Agilent capture. Subsequently, 20 DNA samples from 17 different families, including four 
patients with known mutations were sequenced using lllumina Genome Analyzer llx next-generation-sequencing 
(NGS) platform. Different filtering approaches were applied to identify the genetic defect. The most likely disease 
causing variants were analyzed by Sanger sequencing. Co-segregation and sequencing analysis of control samples 
validated the pathogenicity of the observed variants. 

Results: The phenotype of the patients included retinitis pigmentosa, congenital stationary night blindness, Best 
disease, early-onset cone dystrophy and Stargardt disease. In three of four control samples with known genotypes 
NGS detected the expected mutations. Three known and five novel mutations were identified in NR2E3, PRPF3, EYS, 
PRPF8, CRB1, TRPM1 and CACNA1F. One of the control samples with a known genotype belongs to a family with 
two clinical phenotypes (Best and CSNB), where a novel mutation was identified for CSNB. In six families the 
disease associated mutations were not found, indicating that novel gene defects remain to be identified. 

Conclusions: In summary, this unbiased and time-efficient NGS approach allowed mutation detection in 75% of 
control cases and in 57% of test cases. Furthermore, it has the possibility of associating known gene defects with 
novel phenotypes and mode of inheritance. 

Keywords: NGS, retinal disorders, diagnostic tool. 



Background 

Inherited retinal disorders affect approximately 1 in 2000 
individuals worldwide [1]. Symptoms and associated phe- 
notypes are variable. In some groups the disease can be 
mild and stationary such as in congenital stationary night 
blindness (CSNB) or achromatopsia (ACHM), whereas 
other disorders are progressive leading to severe visual 
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impairment such as in rod-cone dystrophies, also known as 
retinitis pigmentosa (RP) or cone and cone-rod dystro- 
phies. The heterogeneity of these diseases is reflected in the 
number of underlying gene defects. To date more than 150 
genes have been implicated in different forms of retinal dis- 
orders http://www.sph.uth.tmc.edu/Retnet/home.htm and 
yet in a significant proportion of patients the disease caus- 
ing mutation could not be identified, suggesting additional 
novel genes that remain to be discovered. Furthermore, 
recent studies have outlined that distinct phenotypes can 
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be related to the dysfunction of the same gene [2-4]. 
Furthermore, there may be additional phenotype-genotype 
associations that are still not recognized. The state-of-the- 
art phenotypic characterization including precise family 
history and functional as well as structural assessment (i.e. 
routine ophthalmic examination, perimetry, color vision, 
full field and multifocal electroretinography (ERG), fundus 
autofluorescence (FAF) imaging and optical coherence 
tomography (OCT)) allows targeted mutation analysis for 
some disorders. However, in most cases of inherited retinal 
diseases, similar phenotypic features can be due to a large 
number of different gene defects. 

Various methods can be used for the identification of 
the corresponding genetic defect. All these methods have 
advantages and disadvantages. Sanger sequencing is still 
the gold-standard in determining the gene defect, but due 
to the heterogeneity of the disorders it is time consuming 
and expensive to screen all known genes. Mutation detec- 
tion by commercially available APEX genotyping microar- 
rays (ASPER Ophthalmics, Estonia) [5,6] allows the 
detection of only known mutations. In addition, a separate 
microarray has been designed for each inheritance pattern, 
which tends to escalate the costs especially in simplex 
cases, for which inheritance pattern cannot be predeter- 
mined. Indirect methods with single nucleotide poly- 
morphism (SNP) microarrays for linkage and 
homozygosity mapping are also powerful tools, which has 
proven its reliability in identifying novel and known gene 
defects [7-12]. However, in case of homozygosity mapping 
the method can only be applied to consanguineous 
families or inbred populations. To overcome these chal- 
lenges, we designed a custom sequencing array in colla- 
boration with a company (IntegraGen, Evry, France) to 
target all exons and part of flanking sequences for 254 
known and candidate retinal genes. This array was subse- 
quently applied through NGS to a cohort of 20 patients 
from 17 families with different inheritance pattern and 
clinical diagnosis including RP, CSNB, Best disease, early- 
onset cone dystrophy and Stargardt disease. 

Methods 

Clinical investigation 

The study protocol adhered to the tenets of the Declara- 
tion of Helsinki and was approved by the local Ethics 
Committee (CPP, He de France V). Informed written 
consent was obtained from each study participant. Index 
patients underwent full ophthalmic examination as 
described before [13]. Whenever available, blood sam- 
ples from affected and unaffected family members were 
collected for co-segregation analysis. 

Previous molecular genetic analysis 

Total genomic DNA was extracted from peripheral blood 
leucocytes according to manufacturer's recommendations 



(Qiagen, Courtaboeuf, France). DNA samples from some 
patients with a diagnosis of RP were first analyzed and 
excluded for known mutations by applying commercially 
available microarray analysis (arRP and adRP ASPER 
Ophthalmics, Tartu, Estonia). In some cases, pathogenic 
variants in EYS, C2orf71, RHO, PRPF31, PRPH2 and RP1 
were excluded by direct Sanger sequencing of the coding 
exonic and flanking intronic regions of the respective 
genes [13-17]. Conditions used to amplify PRPH2 can be 
provided on request. 

Molecular genetic analysis using NGS 

A custom-made SureSelect oligonucleotide probe library 
was designed to capture the exons of 254 genes for dif- 
ferent retinal disorders and candidate genes according 
to Agilent's recommendations (Table 1). These genes 
include 177 known genes underlying retinal dysfunction 
(http:/ /www.sph.uth.tmc.edu/ retnet/ sum-dis.htm, Octo- 
ber 2010, Table 1) and 77 candidate genes associated 
with existing animal models and expression data (Table 
2). The eArray web-based probe design tool was used 
for this purpose https://earray.chem.agilent.com/earray. 
The following parameters were chosen for probe design: 
120 bp length, 3x probe-tiling frequency, 20 bp overlap 
in restricted regions, which were identified by the imple- 
mentation of eArray's RepeatMasker program. A total of 
27,430 probes, covering 1177 Mb, were designed and 
synthesized by Agilent Technologies (Santa Clara, CA, 
USA). Sequence capture, enrichment, and elution were 
performed according to the manufacturer's instructions 
(SureSelect, Agilent). Briefly, three ug of each genomic 
DNA were fragmented by sonication and purified to 
yield fragments of 150-200 bps. Paired-end adaptor oli- 
gonucleotides from Illumina were ligated on repaired 
DNA fragments, which were then purified and enriched 
by six PCR cycles. 500 ng of the purified libraries were 
hybridized to the SureSelect oligo probe capture library 
for 24 h. After hybridization, washing, and elution, the 
eluted fraction underwent 14 cycles of PCR-amplifica- 
tion. This was followed by purification and quantifica- 
tion by qPCR to obtain sufficient DNA template for 
downstream applications. Each eluted-enriched DNA 
sample was then sequenced on an Illumina GAIIx as 
paired-end 75 bp reads. Image analysis and base calling 
was performed using Illumina Real Time Analysis 
(RTA) Pipeline version 1.10 with default parameters. 
Sequence reads were aligned to the reference human 
genome (UCSC hgl9) using commercially available soft- 
ware (CASAVA1.7, Illumina) and the ELANDv2 align- 
ment algorithm. Sequence variation annotation was 
performed using the IntegraGen in-house pipeline, 
which consisted of gene annotation (RefSeq), detection 
of known polymorphisms (dbSNP 131, 1000 Genome) 
followed by mutation characterization (exonic, intronic, 
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Table 1 Known retinal disease genes Table 1 Known retinal disease genes (Continued) 



Number 


Gene name 


48 


DMD 


1 


ABCA4 


49 


DPP3 


2 


ABCC6 


50 


EFEMPl 


3 


ADAM9 


51 


ELOVL4 


4 


AHI1 


52 


ERCC6 


5 


AIPL1 


53 


EYS 


6 


ALMS1 


54 


CA ft/It £•! A 

FAM161A 


7 


ARL6 


55 


FBLN5 


8 


ARMS2 


56 


FSCN2 


9 


ATXN7 


57 


FZD4 


10 


BBS 10 


58 


GNAT1 


1 1 


BBS 12 


59 


GNAT2 


12 


BBS2 


60 


GPR98 


13 


BBS4 


61 


GRK1 


14 


BBSS 


62 


GRM6 


15 


BBS7 


63 


GUCA1A 


16 


BBS9 


64 


GUCA1B 


17 


BEST1 


65 


GUCY2D 


18 


C1QTNF5 


66 


LJAAfKl 1 

FtmLNl 


19 


C2 


67 


HTRA1 


20 


C2orf71 


68 


IDH3B 


21 


C3 


69 


IMPUH i 


22 


CA4 


70 


IMPG2 


23 


CABP4 


71 


INPP5E 


24 


CACNA1F 


72 


INVS 


25 


CACNA2D4 


73 


IQCB1 


26 


CC2D2A 


74 


JAG1 


27 


CDH23 


75 


KCNJ13 


28 


CDH3 


76 


KCNV2 


29 


CEP290 


77 


KLHL7 


30 


CERKL 


78 


LCA5 


31 


CFB 


79 


LRAT 


32 


CFH 


80 


LRP5 


33 


CHM 


81 


MERTK 


34 


CLN3 


82 


MFRP 


35 


CLRN1 


83 


MKKS 


36 


CNGA1 


84 


MKS1 


37 


CNGA3 


85 


MIND i 


38 


CNGB1 


86 


M/ A/Do 


39 


CNGB3 


87 


MT-AP6 


40 


CNNM4 


88 


MTND2 


41 


COL11A1 


89 


MTND5 


42 


COL2A1 


90 


MTND4 


43 


COL9A1 


91 


MY07A 


44 


CRB1 


92 


NDP 


45 


CRX 


93 


NPHP1 


46 


CYP4V2 


94 


NPHP3 


47 


DFNB31 


95 


NPHP4 
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Table 1 Known retinal disease genes (Continued) 



96 


NR2E3 


97 


NRL 


98 NYX 


99 OAT 


100 


OFD1 


101 


OPA1 


102 


OPA3 


103 


OPN1LW 


104 


OPN1MW 


105 


OPNISw 


106 


OTX2 


107 


PANK2 


1 08 PAX2 


109 


PC DM 5 


110 


PCDH21 


1 1 1 PDE6A 


112 


PDE6B 


113 


PDE6C 


114 


PDE6G 


115 


PDZD7 


116 


PEX1 


117 


PEX2 


1 1 8 PEX7 


1 1 9 PGK1 


120 


PHYH 


121 


PITPNM3 


122 


PRCD 


123 


PROM1 


124 


PRPF3 


125 


PRPF31 


1 26 PRPF8 


127 


PRPH2 


128 


RAX2 


129 


RBI 


130 


RBP3 


131 


RBP4 


132 


RD3 


133 


RDH12 


134 


RDH5 


135 


RGR 


136 


RGS9 


137 


RGS9BP 


138 


RHO 


139 


RIMS1 


140 


RLBP1 


141 


ROM 7 


142 


RP1 
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Table 1 Known retinal disease genes (Continued) 



143 


RP1L1 


144 


RP2 


145 


RP9 


146 


RPE65 


147 


RPGR 


148 


RPGRIP1 


149 


RPGRIP1L 


150 


RSI 


151 


SAG 


152 


SDCCAG8 


153 


SEMA4A 


154 


SLC24A1 


155 


SNRNP200 


156 


SPATA7 


157 


TEAD1 


158 


TIMM8A 


159 


TIMP3 


160 


TLR3 


161 


TLR4 


162 


TMEM126A 


163 


TOPORS 


164 


TREX1 


165 


TRIM32 


166 


TRPM1 


167 


TSPAN12 


168 


TTC8 


1 69 TIP A 


170 


TULP1 


171 


UNC1 19 


172 


USH1C 


173 


USH1G 


174 


USH2A 


175 


VCAN 


176 


WFS1 


177 


ZNF513 



silent, nonsense etc.). For each position, the exomic fre- 
quencies (homozygous and heterozygous) were deter- 
mined from all the exomes already sequenced by 
IntegraGen and the exome results provided by HapMap 
project. 

Investigation of annotated sequencing data 

We received the annotated sequencing data in the form 
of excel tables. On average 946 SNPs and 83 insertions 
and deletions were identified for each sample (Figure 1). 
By using the filtering system, we first investigated var- 
iants (nonsense and missense mutations, intronic 
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Table 2 Candidate genes for retinal disorders 



Number 


Gene name 


Reason 


References 


1 


ADCY1 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


2 


ANKRD33 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


3 


ANXA2 


Promotion of choroidal neovascularization 


[36] 


A 


ARL13B 


Cilia protein, mutations lead to Joubert Syndrome 


[37] 


5 


BMP7 


Regulation of Pax 2 in mouse retina 


[38] 


6 


BSG 




Thierry Leveillard personal commmunication 


7 


CAMK2D 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


8 


CCDC2od 


Moditier tor Bob 


[39,40] 


9 


CLCN7 


Cln7-/- mice severe osteopetrosis and retinal degeneration 


[41] 


1 0 


COL4A3 


Alport syndrome, with eye abnormalities 


[42,43] 


1 1 


COL4A4 


Alport syndrome, with eye abnormalities 


[42,44] 


1 2 


COL4A5 


Alport syndrome, with eye abnormalities 


[42,45] 


13 


CUBN 




Personal communication Renata Kozyraki 


l -1 


CYP1B1 


glaucoma 


[46] 


1 5 


DOHH 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


1 6 


DSC AM L 7 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


1 7 


ESRRB 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


1 8 


FIZ1 


Interactor of NRL 


[47] 


1 9 


GJA9 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


20 


GNAZ 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


21 


GNGT1 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


22 


GPR152 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


23 


HCN1 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


2-1 


LI C A TO it /I 

nc/l TR5A 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


25 


LJICT1 LJ 1 f 

HiSTlHlC 


Expressed in retina 


Expression databases 


26 


IMPG1 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


27 


INSL5 


diff. Expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


28 


KCNB1 


diff. expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


29 


KCTD7 


Expressed in retina 


Expression databases 


30 


LASS4 


diff. expression rd1 mouse 


Chalmel et al., manuscript in preparatiom 


31 


LRIT2 


diff. expression rd1 mouse 


Chalmel et al., manuscript in preparatiom Rd 1 mouse 


32 


LRP2 




Personal communication Renata Kozyraki 


33 


A A A O 1 1 1 1 

MAd2iL1 


diff. expression Rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


34 


MAP2 


diff. expression rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


35 


MAS1 


Degeneration of cones due to expression of Mas1 


[48] 


36 


MAST2 


diff. expression rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


37 


a jinn ^ 

MPP4 


diff. expression rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


38 


MYOC 


glaucoma 


[49] 


39 


NDUrAl2 


diff. expression rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


'10 


NtURUUl 


BETA2/NeuroD1 -/- mouse: photoreceptor degeneration 


rem 


41 


N0S2 


glaucoma 


[51] 


42 


NXNL1 


Rod-derived cone viability factor 


[52] 


43 


NXNL2 


Rod-derived cone viability factor 2 


[53] 


44 


0PN1MW2 


Cone opsin, medium-wave-sensitive2 


[54] 


45 


OPTN 


glaucoma 


[55] 


46 


PFKFB2 


diff. expression rd 1 mouse 


Chalmel et al., manuscript in preparatiom 


47 


PIAS3 


Rod photoreceptor development 


[56] 
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Table 2 Candidate genes for retinal disorders (Continued) 



48 


PKD2L1 


Diff. expression in human retinal detachment 


Delyfer et al. 201 1 submitted 


49 


PLEKHA1 


Age-related macular degeneratiom 


[57] 


50 


PPEF2 


diff. expression rdl mouse 


Chalmel et al., manuscript in preparatiom 


51 


RAB8A 


Interacts with RPGR, role in cilia biogenesis and maintenance 


[58] 


52 


RABGEF1 


diff. expression rdl mouse 


Chalmel et al., manuscript in preparatiom 


53 


RCVRN 


diff. expression rdl mouse 


Chalmel et al., manuscript in preparatiom 


54 


RGS20 


diff. expression rdl mouse 


Chalmel et al., manuscript in preparatiom 


55 


RNF144B 


diff. expression rdl mouse 


Chalmel et al., manuscript in preparatiom 


56 


RORB 


Rod photoreceptor development in mice 


[59] 


57 


RXRG 


Retinoic acid receptor, highly expressed in the eye 


Expression databases 


58 


SGI PI 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


59 


SLC16A8 


Altered visual function in ko-mice 


[60] 


60 


SLC17A7 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


61 


STAM2 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


62 


STK35 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


63 


STX3 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


64 


SV2B 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


65 


TBC1D24 


diff. expression rd 1 mouse 


Chalmel et al, manuscript in preparatiom 


66 


THRB 


Essential for M-cone development in rodents 


[61] 


67 


TMEM216 


Cilia protein, mutations lead to Joubert and Meckel syndrome 


[62] 


68 


TMEM67 


Cilia protein, mutations lead to Joubert 


[63] 


69 


TRPC1 


diff. expression rd 1 mouse 


diff. expression Rd 1 mouse 


70 


UHMK1 


diff. expression rdl mouse 


diff. expression Rdl mouse 


71 


VSX1 


Stimulator for promoter NXNL1 


[64] 


72 


VSX2 


Stimulator for promoter NXNL1 


[64] 


73 


WDR17 


diff. expression rdl mouse 


diff. expression Rdl mouse 


74 


WDR31 


diff. expression Nxnh-/- mouse 


[65] 


75 


WISP1 


diff. expression rdl mouse 


Chalmel et al, manuscript in preparatiom 


76 


XIAP 


Protects photoreceptors in animal models of RP 


[66] 


77 


ZDHHC2 


diff. expression Rdl mouse 


Chalmel et al, manuscript in preparatiom 



variants located +/- 5 apart from exon), which were 
absent in dbSNP and NCBI databases http://ncbi.nlm. 
nih.gov/. In the absence of known gene defects or puta- 
tive pathogenic variants (see below) in the first step, we 
selected known genes, which were previously clinically 
associated including variants present in dbSNP and 
NCBI databases (Figure 1). Each predicted pathogenic 
variant was confirmed by Sanger sequencing. 

Assessment of the pathogenicity of variants 

Following criteria were applied to evaluate the patho- 
genic nature of novel variations identified by NGS: 1) 
stop/frameshift variants were considered as most likely 
to be disease causing; 2) co-segregation in the family; 3) 
absence in control samples; 4) for missense mutations 
amino acid conservation was studied in the UCSC Gen- 
ome Browser http://genome.ucsc.edu/ across species 
from all different evolutionary branches. If the amino 



acid residue did not change it was considered as "highly 
conserved". If a different change was seen in fewer than 
five species and not in the primates then it was consid- 
ered as "moderately conserved" and if a change was pre- 
sent in 5-7, it was considered as "weakly conserved", 
otherwise the amino acid residue was considered as "not 
conserved", 5) pathogenicity predictions with bioinfor- 
matic tools (Polyphen: Po/ymorphism Phenotyping, 
http://genetics.bwh.harvard.edu/pph/ and SIFT: Sorting 
Intolerant From Tolerant, http://blocks.fhcrc.org/sift/ 
SIFT.html) if at least one of the program predicted the 
variant to be possibly damaging, it was considered to be 
pathogenic; 6) presence of the second mutant allele in 
the case of autosomal recessive inheritance. Mutations 
were described according to the HGVS website http:// 
www.hgvs.org/mutnomen. In accordance with this 
nomenclature, nucleotide numbering reflects cDNA 
numbering with +1 corresponding to the A of the ATG 
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IntegraGen 
Variant detection: 



946 SNPsand83 InDels 



•dbSNP and NCBI filtering 
•Including only exonic and +/-5 intronic variants 
•Including only missense and nonsense mutations 
•Exchange never present in 1000 genomes: 



11 variants 

4 

•Known gene defect? 

•mode of inheritance 

•protein function 

'polyphen, sift: ) 

■ 

1-5 variants 




•No dbSNP and NCBI filtering 
•For ar cases: Exchanges can be 
het present in 1000 genomes 
•Known gene defect? 
i 

1-5 variants 



•Conformation by Sanger sequencing 
•Co-segregation analysis 
•Exclusion in control samples 



Figure 1 Flow chart of variant analysis. IntegraGen provided the results in form of excel tables. For each sample on average 946 SNPs and 83 
inDels were detected, of which 11 represent missense, nonsense or putative splice site mutations, which were absent in dbSNB, NCBI and 1000 
genome databases. Of those 1-5 variants were predicted to be pathogenic. In case where none of the variants were predicted to be pathogenic, 
dbSNB, NCBI and 1000 genome databases were included to detect mutations referenced with an rs-number. Co-segregation analysis was 
performed in families with putative pathogenic variants. 



translation initiation codon in the reference sequence. 
The initiation codon is codon 1. The correct nomencla- 
ture for mutation was checked applying Mutalyzer 
http:/ /www.lovd.nl/ mutalyzer/. 

Results 

The overall sequencing coverage of the captured regions 
was 98.4% and 90.4% for a lx and a lOx coverage respec- 
tively. The overall sequencing depth was > 120x. The 
number of reference and variant sequences detected by 
NGS, reflected the correct zygosity state of the variant; 
on average if 50% of the sequences represented the var- 
iant, then a heterozygous state was called, while if 100% 
of the sequences represented the variant, then a homozy- 
gous or hemizygous state was annotated by IntegraGen. 



Validation of the novel genetic testing tool for retinal 
disorders 

To validate the novel genetic testing tool for retinal dis- 
orders, we used four DNA samples from families, in 
which we had previously identified different types of 
mutations by Sanger sequencing: one 1 bp duplication 
and one 1 bp deletion in PRPF31 and missense muta- 
tions in TRPM1 and BEST1 (Table 3). Three of the four 
mutations were detectable by NGS, whereas the deletion 
in PRPF31 was not identified. To validate if this was due 
to a technical problem of deletion detection in general 
or low coverage at this position, the sequencing depth 
was investigated in detail. Indeed the coverage at this 
position reflected by the mean depth was only ~l-6 for 
all samples. This indicates that although the coverage in 
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Table 3 Patients with known mutations used to validate the novel genetic approach for retinal disorders 



Index 


Phenotype 


Gene 


Mutation 


Allele 


Read reference 

IN vj j 


Read variant 


Mutation detected by 


Mean 
depth 


CIC00034, 

F28 


adRP 


PRPF31 


c.666dup 
p 

l223YfsX56 


het 


11 


13 


yes 


21.3-22.5 


CIC00140, 
F108 


adRP 


PRPF31 


c.997delG 
P- 

E333SfsX5 


het 






no 


5.0-5.2 


CIC00238, 
F165 


arCSNB 


TRPM1 


C.1418G > 
C 

p.R473P 


homo 


0 


38 


yes 


36.7 


CIC00707, 
F470 


Best and adCSNB see 
Table 5 


BEST1 


C.73C > T 
p.R25W 


het 


40 


38 


yes 


99.4 



general was very good, specific probes used here need to 
be redesigned to improve the capture for specific exons. 

Detection of known and novel mutations 

Some of the patients from the 14 families with no 
known gene defect were previously excluded for known 
mutations using microarray analysis and by Sanger 
sequencing in the known genes EYS, C2orf71, RHO, 
PRPF31, PRPH2 and RP1. Other samples were never 
genetically investigated. In four DNA samples known 
mutations were detected (Table 4) from three different 
families with autosomal dominant (ad) or recessive (ar) 
RP. All mutations co-segregated with the phenotype 
(Figure 2). In seven samples, novel mutations in known 
genes were identified. These mutations co-segregated 
with the phenotype from five different families with 
adCSNB, x-linked incomplete CSNB, adRP, arRP and x- 
linked RP (Table 5, Figures 3 and 4). One of the cases 
from these five families was also used as a control for 
Best disease carrying a known BEST1 mutation (Table 
3). In addition to the Best phenotype, ERG-responses of 
this patient resembled those of complete CSNB, i.e. 
showing selective ON-bipolar pathway dysfunction. This 
phenotype was independent of the Best phenotype (Fig- 
ure 3). The most likely disease causing mutation 



detected by NGS was a novel heterozygous TRPM1 
mutation (Table 4, Figure 3). 

Unsolved cases 

In six of the 14 families with Stargardt disease, adRP, 
adCD with postreceptoral defects, arRP, early onset 
arCD with macrocephaly and mental retardation 
described in affected sister and x-linked cCSNB, the dis- 
ease associated mutations remain to be elucidated or 
validated (Table 6, Figure 5). 

Discussion 

By using NGS in 254 known and candidate genes we 
were able to detect known and novel mutations in 57% 
of families tested. In order to achieve this goal, we 
applied a rigorous protocol (Figure 1). To our knowl- 
edge, this is the first report using NGS to investigate all 
inherited retinal disorders at once. In a study restricted 
to adRP, Bowne and co-workers used a similar approach 
including 46 known and candidate genes for adRP [18]. 
All their cases had previously been screened and 
excluded for most of the known genes underlying adRP. 
The authors were able to identify known or novel muta- 
tions in five out of 21 cases in genes not included in a 
pre-screening [18]. This added five patients to their 



Table 4 Detection of known mutations by using the novel genetic approach for retinal disorders 



Index 


Phenotype 


Pre-screening 


Gene 


Mutation 


Allele 
State 


Read 

reference 

NGS 


Read 

variant 

NGS 


Reference 


Mutation verified by 
Sanger and co- 
segregation 


CIC00019, 
F16 


adRP 


Linkage, RHO, 
PRPF31, PRPH2, 
RP] 


PRPF3 


C.1481C > T 
p.T494M 


het 


25 


22 


[67] 


yes 


CIC0000893, 
F574 


adRP 


RHO, PRPF31, 
PRPH2, RP1 


NR2E3 


C.166G > A 
p.G56R 


het 


5 


3 


[68] 


yes 


CIC000128, 
F100 


arRP, 
consang. 




EYS 


c.408_423del 
p.N137VfsX24 


homo 




179 


[1 3,69] 


yes 


CIC0000943, 
F100 


arRP, 
consang 




EYS 


c.408_423del p. 
N137VfsX24 


homo 


0 


193 


[1 3,69] 


yes 
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Figure 2 Detection of known mutations by NGS in 254 retinal genes. The index patient 19 of family 16 with adRP revealed the p.T494M 
mutations in PRPF3, which co-segregates with the phenotype. Two family members never clinically investigated from the last generation (984 
and 1 167 carrying a question mark) were reported to be not affected but carried the mutation. They may develop the phenotype at a later 
stage. In addition variability of the phenotype of this mutation was documented [35], Two patients, 128 and 943 of family 100 with arRP from 
Jewish origin revealed the known EYS mutation p.N137VfsX24, which was found in all screened affected family members. The index patient 893 
of family 574 showed the previously described NR2E3 p.G56R mutation, which co-segregated with the phenotype. 



adRP cohort with known gene defects, indicating that 
64% of their patients show known mutations with new 
genes still to be discovered in the remaining 36%. The 
current study provides a more exhaustive tool, since it 
incorporates screening of 254 genes implicated in var- 
ious retinal disorders of different inheritance patterns 
and additional candidate genes for these phenotypes. 
With this approach a cohort of both pre-screened and 
unscreened samples, was investigated. The mutation 
detection rate of 57% is high and was never obtained 
before by high throughput screening methods. Further- 
more, this approach is probably less time consuming 
and expensive than existing methods such as direct 
sequencing of all known genes or microarray analysis. 
Of note however is one of the variants detected with the 
NGS approach (i.e. p.V973L exchange in GUCY2D), 



which was not confirmed by direct Sanger sequencing, 
suggesting the possibility of false positive using the high 
throughput screening. Verification by direct Sanger 
sequencing of most likely pathogenic variants is there- 
fore essential to validate NGS data, although the false 
positive rate is assumed to be low (in our study 1/28 
verified sequence variants represented a false positive). 

Overall, the study of 20 subjects from 17 families by 
NGS showed that most of the targeted regions are well 
covered (more than 98%). However, some of the regions 
showed a lower coverage (GC-rich regions) or were not 
captured (repetitive regions). This was for instance the 
case for two genes underlying cCSNB, (i.e. NYX and 
GRM6) and the repetitive region of ORF15 of RPGR. 
For GC-rich regions the capture design could be 
improved in the future by modifying NGS chemistry, as 
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Table 5 Detection of novel mutations by using the novel genetic approach for retinal disorders 

Index Phenotype Pre- Gene Mutation Allele Read Read Mutation Conservation Polyphen Sift 

screening State reference variant verified by 

NGS NGS Sanger and 



co- 
segregation 



CIC00707, 
F470 


adCSNB 
and Best 
see Table 3 


RHO, 

PDE6B, 

GNATl 


TRPM1 


C.1961A 

> C 

p.H654P 


net 


39 


38 


yes 


moderately 
conserved 


possibly 
damaging 


tolerated 


CIC000348, 
F232 


adRP, mild 


RHO, 
PRPF3 1, 
PRPH2, 
RPl, adRP 
chip 


PRPF8 


C.6992A 
> G 
P- 

E2331G 


het 


13 


10 


yes 


moderately 
conserved 


possibly 
damaging 


affect 

protein 

function 


CIC000346, 
F232 


adRP 




PRPF8 


C.6992A 
> G 

p.E2331G 


het 


5 


9 


yes 


moderately 
conserved 


possibly 
damaging 


affect 

protein 

function 


CIC000347, 
F232 


as 

adRP 




PRPF8 


C.6992A 

> G 

p.E2331G 


het 


15 


17 


yes 


moderately 
conserved 


possibly 
damaging 


affect 

protein 

function 


CIC04240, 
F2025 


arRP, 
consang., 
detailed 
clinic in 
[70] 


RSI 


CRB1 


C.2219C 
> T 

p.S740F 


homo 


2 


194 


yes 


highly 
conserved 


probably 
damaging 


affect 

protein 

function 


CIC00199, 
F146 


adRP or x- 
linked RP 
with 
affected 
carrier 


RHO, 
PRPF31, 
PRPH2, 
RPl, adRP 
chip 


RPGR 


C.248-2A 
> G 
splice 
defect 


hetero 


30 


22 


yes 


conserved 
splice site 


n.a. 


n.a. 


CIC04094, 


icCSNB 




CACNA1F 


C.973C > 


hemi 


0 


28 


yes 


n.a. 


n.a. 


n.a. 



F1915 T 

p.Q325X 



it was successfully achieved for Sanger sequencing using 
different additives, which improved the amplification 
and subsequent sequencing. If repetitive regions like 
ORF15 of RPGR remain problematic for sequencing by 
NGS, direct Sanger sequencing of these targets might be 
the first screening of choice; in particular for disorders 
caused only by a few gene defects such as CSNB, and 
xl-RP. 

By applying NGS sequencing to our retinal panel, 
known and novel mutations were detected in different 
patients. We believe that our diagnostic tool is particu- 
larly important for heterogeneous disorders like RP, for 
which many gene defects with different prevalence have 
been associated to one phenotype. It also allows the 
rapid detection of novel mutations in minor genes 
which are often not screened as a priority by direct San- 
ger sequencing. This was the case in our study for three 
individuals from one family with adRP in which NGS 
detected a novel PRPF8 mutation in both affected and 
one unaffected family member (Table 4, Figure 4). In 
this family, the RP phenotype is mild and therefore it is 
possible that the unaffected member may develop symp- 
toms later in life or alternatively it may be a case of 
incomplete penetrance as reported for another splicing 



factor gene, PRPF31 and recently for PRPF8 as well 
[19-22]. Interestingly, a novel TRPM1 mutation was 
identified in a patient with adCSNB, a gene previously 
only associated with arCSNB [23-26]. This is the first 
report of a TRPM1 mutation co-segregating with ad 
Schubert-Bornschein type complete CSNB. Since the 
location of this mutation is not different compared to 
other mutations leading to arCSNB, it is not quite clear 
how TRPM1 mutations might lead to either ad or 
arCSNB. Functional investigations are needed to validate 
the pathogenicity of this variant. Furthermore, this find- 
ing suggests that TRPM1 heterozygous mutation carriers 
from arCSNB families should be investigated by electro- 
retinography to determine whether they display similar 
retinal dysfunction as in affected members of the pre- 
sented adCSNB family. Detection of a novel RPGR splice 
site mutation in family 146 presented a challenge. The 
actual disease causing change was concealed under a 
wrongly annotated rs62638633, which had previously 
been clinically associated to RP by a German group 
http://www.ncbi.nlm.nih.gov/sites/varvu7gen- 
e=6103&rs=62638633, (personal communication, Mar- 
kus Preising). These observations indicate that the 
stringent filtering we applied initially can mask those 
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Fam 470: BEST1: Ml: c.73C>T p.R25W; TRPM1: M2:c.l961A>C p.H654P 
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Figure 3 Best disease and CSNB co-segregating in one family, a) Sanger and NGS detected in all patients with Best disease a BEST1 
mutation, b) NGS detected in all patients with a cCSNB phenotype a novel TRPM1 mutation, c) Fundus colour photographs (above) and fundus 
autofluorescence (below) of patient 707 showing multiple yellow deposits within the posterior pole which are hyper autofluorescent d) Electro- 
oculogram of patient 707 showing no slight rise after illumination in keeping with the diagnosis of Best disease e) Full Field Electroretinogram of 
patient 707 showing ON-bipolar cell pathway dysfunction in keeping with the diagnosis of cCSNB. 



referenced disease causing variants. Bearing this in mind 
one can still first investigate unknown variants, but 
should then examine dbSNP for referenced variants 
either described to be disease causing, having a low 
minor allele frequency or present in interesting candi- 
date genes. An accurate discrimination of non-patho- 
genic polymorphisms versus disease causing 
polymorphism in SNP databases is warranted to resolve 
this challenge. 

In six families from the investigated cohort the disease 
causing mutations still remain to be identified. In the 
Stargardt patient with no pathogenic ABCA4 mutations 
two variants in CFH were detected, one of which 
(rsl061170) had previously been reported to predispose 
to age related macular degeneration (AMD) [27-29]. The 
second CFH change is a novel variant, affecting a highly 
conserved residue, not found in NGS data from the other 
19 samples and never associated with a disease. The 



variants co-segregated in the only available family mem- 
bers, which were the patient's parents. Apart from the 
association with AMD, CFH mutations have been pre- 
viously associated with renal diseases, the most common 
being membranoproliferative glomerulonephritis and 
hemolytic uremic syndrome, which can be also associated 
with an eye phenotype [30,31]. No renal dysfunction was 
present in our patient. To validate if the two variants 
identified in CFH are indeed disease causing, the DNA 
samples from other available family members for co-seg- 
regation analysis as well as characterization of functional 
consequences of the novel variant are needed. One 
patient with complete CSNB had an affected nephew and 
thus x-linked inheritance was assumed. However, neither 
Sanger nor NGS detected a mutation in the only known 
x-linked gene, NYX, causing cCSNB. To exclude reces- 
sive inheritance TRPM1 and GRM6 were investigated in 
detail. Indeed the patient carried a novel heterozygous 
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Figure 4 Detection of novel mutations using NGS in 254 retinal genes. Novel mutations in PRPF8, CRB1, RPGR and CACNA1F co-segregated 
in affected and asymptomatic carriers with the adRP, arRP, x-linked dominant and X-liked icCSNB phenotypes respectively. Asymptomatic 
individuals are marked with a question mark. 



TRPM1 variant, which affects a highly conserved amino 
acid and was not identified in the other 19 samples inves- 
tigated here (Table 6). However, direct Sanger sequen- 
cing of lower covered regions did not identify a second 
mutation in this gene. Similarly no mutations in GRM6 
were identified. These findings outline the need for addi- 
tional family members to determine, through co-segrega- 
tion, the pathogenicity of the numerous variants 
identified by NGS. This was also true for two other 
families with nonsense mutations in CUBN (Fam795) 
and RP1L1 (Fam761) (Table 6). The nonsense mutation 
in CUBN, co-segregated with the phenotype in most of 
the family members (Figure 5). Had we not had access to 
additional family members, we might have retained this 
gene defect as the underlying cause for adCD and consid- 
ered CUBN as a new gene involved in adCD. None of the 
other putatively pathogenic mutations identified in 
CUBN, TRPM1 and GUCY2D co-segregated with the 
phenotype in this family (Table 6, Figure 5). RP1L1 was 
already a candidate for adRP [32] but was previously 



associated with occult macular dystrophy [33]. In our 
study, this variant did not co-segregate with the pheno- 
type in other affected family members (data not shown). 

This NGS study ended with six genetically unresolved 
families, which can be further investigated with whole 
exome sequencing. Although, no clear information 
about the actual percentage of missing gene defects 
underlying each group of inherited retinal disorders 
exists, previous studies have reported that in many cases 
the genetic cause still needs to be determined [18,34]. 
Whole exome sequencing approaches allow the detec- 
tion of both, novel and known gene defects, but also 
generate numerous variants and therefore require the 
inclusion of more than one DNA sample for each family 
to rapidly exclude non-pathogenic variants. Due to the 
higher costs of exome sequencing for one sample com- 
pared to targeted sequencing, we propose to initially 
perform targeted sequencing in the index patient and 
proceed only after exclusion of a known gene defect to 
whole exome sequencing. 
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Table 6 Patients with unsolved genotype and unlikely disease causing mutations 
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Table 6 Patients with unsolved genotype and unlikely disease causing mutations (Continued) 
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Figure 5 Detection of novel mutation by using NGS in 254 retinal genes. Family 795 reveals autosomal dominant cone dystrophy with 
post-receptoral defects. Four putative disease causing mutations were investigated on the basis of co-segregation. However, none of them co- 
segregated in all affected family members with the phenotype and thus are not considered to be disease causing. Individuals marked with a star 
were clinically investigated, patients with a question mark are asymptomatic and patients with a plus sign show high myopia. 



Conclusions 

In summary, our diagnostic tool is an unbiased time 
efficient method, which not only allows detecting known 
and novel mutations in known genes but also potentially 
associates known gene defects with novel phenotypes. 
This genetic testing tool can now be applied to large 
cohorts of inherited retinal disorders and should rapidly 



deliver the prevalence of known genes and the percen- 
tage of cases with missing genetic defect for underlying 
forms of retinal disorders. 



List of abbreviations 

ad: autosomal dominant; ar: autosomal recessive; as: asymptomatic; het: 
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